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A method which is well suited for implementation on a digital computer 
is presented for the solutions o] active circuits. Unlike many state variable 
approaches the state vector is defined as the set oj voltages which exist between 
certain nodes and the reference node. An advantage of this approach is that 
degeneration in the order oj complexity oj the network caused by capacitance 
loops is handled automatically. Any type oj controlled source can be specified. 
From the basic algorithm the circuit is specified in matrix form by inspection 
using standard nodal methods, and the solution is obtained by a systematic 
reduction oj this one matrix equation. An upper bound on the order oj 
complexity oj the network is evident from the network topology or Hie 
partitioned form of the original matrix. Inductors are included in this 
approach by considering the equivalent gyrator-capacitor combination. 



I. INTRODUCTION 

State variable techniques presently being used to analyze networks 
require a detailed knowledge of graph theory. 1 - 7 Another method of 
state variable analysis that is based partly on a nodal approach and 
does not require a detailed knowledge of graph theory is very re- 
strictive. 8 The method presented here performs a nodal analysis on a 
transformation of the network in which all magnetic storage elements 
have been replaced by gyrator-capacitor equivalents, and nothing 
more than a basic knowledge of graph theory nomenclature is re- 
quired. The rclmst* network can be transformed to an equivalent 



♦Resistor, capacitor, inductor, mutual inductor, source and ideal transformer. 
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Fig. 1 — Inductor and transformer equivalents. 

resistance capacitance source network using the gyrator-capacitor 
equivalents shown in Fig. 1. Each gyrator shown in Fig. 1 has the 
indefinite admittance parameters 

Ia 

h 
Jc. 

choosing this type of gyrator enables the capacitor value in farads of 
the equivalent pair to be equal to the inductor value in henries. 

Let the number of nodes of a transformed network be n. Using 
Kirchoff's current law, it can be shown that for an n-node RCS net- 
work 

C'V = I - GV (1) 

where I is an (n — l)th ordered column vector representing the currents 
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injected into the nodes, V is an (n — l)th ordered column vector repre- 
senting the voltages between the nodes of the network and the reference 
node. If the transformed network contains I capacitors then the matrix 
C is an (n — l)th ordered symmetric matrix which contains imbedded 
within it I second order indefinite matrices, each having the dimensions 
of farads. Similarly G represents the resistors and has the dimensions of 
mhos, but G may be asymmetric. Node n is the common or ground node 
of the network; for convenience this node is always assumed to have 
capacitors connected to it. 

The objective is to find an upper bound on the rank of the ca- 
pacitance matrix C by partitioning C as described in Section II, and 
reducing the matrix equation (1) containing the partitioned matrix 
C to the rank of C ; this reduction is symbolic and does not take into 
account degenerate cases which can occur. It is shown in Appendix B 
that for all conditions, for any type of circuit, an upper bound on the 
order of complexity of the network (rank of C) can be found from 
the network topology. 

II. PARTITIONING OF THE CAPACITANCE MATRIX 

There are basically four types of voltage source (vs) , the independ- 
ent vs (ivs) , the voltage dependent vs (vdvs) , the current dependent 
vs where the current is through a resistor (cdvsr), and the current 
dependent vs where the current is through a capacitor (cdvsc). It will 
be shown that the only current source (cs) which can effect the parti- 
tioning is the current dependent cs where the current is through a 
capacitor (cdcsc). As a result, any type of cs will be termed simply 
a cs, unless it is a cdcsc. 

The method of partitioning makes the reduction of the matrix 
equation (1) to its rank a simple process. Generally only the voltage 
at a node to which a capacitor is connected can be a state variable 
node. However it is possible to choose a node to which a cdcsc or 
cdvsc is connected as a state variable node instead of one of the 
nodes of the capacitor whose current supplies the dependence, but this 
possibility is avoided automatically in the partitioning method pre- 
sented here. 

The presence of inductors and time-invariant, independent cs's 
forming a cut-set in the original untransformed network causes a 
linear dependence problem in the transformed network. In the trans- 
formed network such a cut-set appears as a capacitor tree with gyra- 
tors only connected to the end nodes of the tree as shown in Fig. 1, 
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and gyrators and perhaps time-invariant independent c.s.'s connected 
to the central node (the gcnode) ; the nodes of this capacitor tree will 
be called the gcset nodes. 

The capacitors in the transformed network can be divided into two 
classes, those connected to the reference node directly or through a 
vs-capacitor chain (the fixed capacitors), and those not so connected 
(the floating capacitors). The ra floating capacitor subgraphs are 
defined as the m unconnected subgraphs obtained from the floating 
capacitor plus imbedded vs graph of the transformed network. 

The partitioning of the capacitance matrix will be related to the 
example of Appendix A in the discussion that follows. Partition the 
matrix C as 
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where 

(i) The nodes nl are all the nodes to which capacitors are connected 
omitting the following nodes: 

(a) A node for each vs imbedded in a capacitor chain (these nodes are 
in the n2 section), but each capacitor must be specified by at least one 
node. 

(b) A node for each of the m floating capacitor subgraphs (these nodes 
being in the n3 section) . 

(c) A node for each gcset which is specified in section n2. 
In the example in Appendix A, nl contains nodes 1 — > 9. 

(ii) The nodes n2 represent: 

(a) A node for each dvs imbedded in a capacitor chain. 

(b) A node free of capacitors for each cdvsc and cdcsc free of capaci- 
tors on at least one node. 
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(c) A node for each gcsbt, 

In the Appendix A example n2 contains nodes 10 and 11. 

(Hi) Section n3 contains a node for each of the m floating capacitor 
subgraphs. 

In the example n3 contains nodes 12 and 13. 

(iv) Section n4 contains the nodes to which only resistors and cs's (but 
not cdcsc) are connected, including a node for each ivs, vdvs or cdcsc 
free of capacitors, other cdvsc's or cdcsc's on both nodes. (The other 
nodes of these sources are specified in section n5). 
In the example n4 has no entries. 

(i 1 ) Section n5 contains all the remaining nodes. These are: 

(a) A node for each ivs. 

(6) A node free of capacitors for each vdvs or cdvsr free of capacitors 
or cdvsc or cdcsc on at least one node. 

In the example n5 contains node 14. 
The rank of the C matrix is nl, and nl = 9 for the example of Appendix A. 
Notice that the presence of capacitance loops in no way alters the method 
of partitioning. 

III. REDUCTION of the circuit description to a minimal form 

Theorem: An upper bound on the order of complexity of a network is the 
order of nl. 

This theorem is proved in Appendix B, where it is shown that every 
row in sections n2, n3, n4, and n5 is linearly dependent on rows in section 
nl; the subspace spanned by sections n2, n3, n4, and n5 is contained in nl. 

The systematic reduction of equation (1) is accomplished by first 
eliminating section n5 by applying the voltage restrictions caused by the 
vs's in section n5. Secondly, section n4 is eliminated using the fact that 
these nodes are free of capacitors. Next, section n3 is eliminated to 
correct the over specification of the floating capacitor subgraphs. 
Finally, the remaining dependencies of the system are caused by the 
dvs's imbedded in capacitive chains, the cdvsc and cdcsc free of capaci- 
tors on at least one node, and a node for each capacitive tree in which a 
gcset has occurred; these dependencies are eliminated with section n2, 
yielding equation (9) of Appendix B. 

Equation (9) of Appendix B can be written as 
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vl = B - Avl (2) 

where 

""vl 

v2 
v = v3 
v4 
_v5. 

and vl is the voltage vector for nodes nl, n2, n,3, etc. 

An example of a solution based on the problem set by Pottle is 
given in Appendix C. 

IV. CONCLUDING REMARKS 

A state variable technique has been described that offers two ad- 
vantages over traditional methods: 

(i) The network can be specified completely by inspection using 
well known nodal techniques with little skill required, the problem 
then becoming one of simple matrix reduction (easily programmed 
for a digital computer). 

(ii) Capacitor loops present no problem and are not even recognized 
as such since the partitioning and matrix reduction are unaltered if 
there are any capacitor loops present. 

The main disadvantages are that currents must always be ex- 
pressed as functions of node voltages and inductors must be replaced 
by gyrators and capacitors; inductor cut-sets must be recognized and 
the circuit redrawn before inductors are eliminated so that the cut- 
set encircles one node only, and this is sometimes inconvenient. 

APPENDIX A 

Example of Partitioning 

For the example of Fig. 2(a), the transformed circuit without in- 
ductors is given in Fig. 2(b). (This is a theoretical problem and the 
circuit has no practical value.) This circuit is described by the equa- 
tions 
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Fig. 2 — Circuit to demonstrate transformation and partitioning. 
where I\ and 7 2 are the unbalance currents due to the vs's and 

Notice that except for degenerate cases (for example, if C 6 = 0) , the 
order of complexity of this network is 9. 



APPENDIX B 

Matrix Reduction 

Consider the partitioned form of equation (1). The section n5 can be 
eliminated as follows: for an ivs of a volts connected between nodes k 
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and I (node k is in section n5, node I is not) 

v k = V t + a. 

For a vdvs or cdvsr connected between nodes k and I, where the vs is 
dependent on the voltage vector v m (each voltage of v m is not in section 
n5). 

v k = Vi + /3v m . 
Thus the system can be reduced to 
nl n2 n3 n4 
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(3) 



Nodes n4 can be eliminated by first writing part of equation (3) as 

"vl 
14 - [G*, G 42 G i3 ] 



v4 = G7A 
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The matrix 



di Gi 2 G13 
G21 G22 G23 

.G31 G32 t/33_ 

has order nl + n2 + n3 and rank no greater than nl + n2. The n3 
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linearly dependent rows and columns can be deleted from equation (5) 
by adding selected rows in the section nl and n2 to rows in the range n3. 
The selection is made as follows: starting with any row in the section n3, 
examine the first entry. If it is nonzero add row 1 to this row. Continue 
along the row, repeating if necessary, until all the entries are zero. 
Proceed for the other dependent rows. Equation (5) can then be written 
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and 



v3 = GV^ktS' - [Gl' n GV»] 
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(7) 



It is a simple process for the reader to prove to himself that elimi- 
nating a node of a floating capacitor subgraph which is part of a 
gcset as described above yields the same result as equating the alge- 
braic sum of the voltages across the capacitors in the gcset to zero 
(analogous to the algebraic sum of the currents entering the inductor 
cut-set node through the inductors adding up to zero) . 

Substituting equation (7) and its derivative into equation (6) we 
obtain 
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(8) 



The total number of restrictions have not yet been placed on the 
network. 

(i) For a dvs imbedded in a capacitor chain or a cdvsc free of capaci- 
tors on one node connected between nodes fc and I, where node k is 
specified in section n2, 

V k = Vt + TV,- 



or 



v k = Vi + i'V,- 

where v,- is the set of voltages upon which the source is dependent. A 
particular voltage of v, may be in any section nl, n2, n3, n4, or n5. 
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(ii) For a cdcsc connected between nodes k and I, the currents 
I k and 7; injected into nodes k and I with the cdcsc removed must be 
modified to 



and 



h + Wi 



h - yv» 



respectively, where 17 has the dimensions of farads. 

(iii) For a gcset with node j of the capacitor tree containing the 
gcset specified in section n2, node j is eliminated as follows: examine the 
entries of row ; of the remaining capacitance matrix. If entry C it i 9* 0, 
subtract Cj,i/Ci, t times row I from row ;, where I = 1, p; p is the order of 
nl + n2. Thus row 3 is reduced to a row of zeros. 

The system can now be written as 

[C]vl = iFl - [G]vl. (9) 

Barring degeneracy, matrix C is nonsingular with rank nl. 

APPENDIX C 

Example of the Method 

For the circuit of Fig. 3 (the example of C. Pottle 9 ), nodes 1, 2, and 3 
are placed in the nl section, and node 4 is placed in the n5 section. Thus, 
by inspection 
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where 



Vi = E; v, = E. 
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Fig. 3 — Example of Appendix C. 

The derivative of the source E must be considered if a capacitor is 
connected to both of its nodes. Clearing out the voltage terms in the 
current array, 



c, -c 4 








c 5 








-c 4 









c 4 ~ 


". 1 
v, 








V-2 


c 3 





h 





cj 


.%. 



G 3 - G 2 
-Go 

-G 3 
. 



G, 

Oi + G 2 



-G x 



-c 3 


G 3 + G\ 








»i 


-<?, 


v 2 





v 3 


G. _ 


-«4_ 



Eliminating y 4 and £> 4 as described in Appendix B, 
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This is as far as we can go symbolically and as far as the method 
takes us. Normally all that remains is a simple inversion of the re- 
maining capacitance matrix, but Pottle chose C\ = C 4 . This makes 
the capacitance matrix singular and so another node must be elimi- 
nated. Eliminating node 1, 
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G~2 G 2 G 3 



G 1 + G 2 + 
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G 3 — G'z 



G 3 + G 4 



G* — G 2 



The vector 



ft] 

can now be expressed explicitly. 
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